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Abstract 

The evolution of the low-lying El strength in proton-rich nuclei is analyzed in the framework of 
the self-consistent relativistic Hartree-Bogoliubov (RHB) model and the relativistic quasiparticle 
random-phase approximation (RQRPA). Model calculations are performed for a series of N=20 
isotones and Z=18 isotopes. For nuclei close to the proton drip-line, the occurrence of pronounced 
dipole peaks is predicted in the low-energy region below 10 MeV excitation energy. From the 
analysis of the proton and neutron transition densities and the structure of the RQRPA amplitudes, 
it is shown that these states correspond to the proton pygmy dipole resonance. 

PACS numbers: 21.10.Gv, 21.30.Fe, 21.60.Jz, 24.30.Cz 
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A number of experimental and theoretical studies of the multipole response of exotic 
nuclei far from stability have been reported in recent years. On the neutron-rich side, in 
particular, the possible occurrence of the pygmy dipole resonance (PDR), i.e. the resonant 
oscillation of the weakly-bound neutron skin against the isospin saturated proton-neutron 
core, has been investigated. The onset of low-lying El strength has been observed not only 
in exotic nuclei with a large neutron excess, e.g. for neutron-rich oxygen isotopes [lj], but also 
in stable nuclei with moderate proton-neutron asymmetry, like 44,48 Ca and 208 Pb ODD 
On the theory side, various models have been employed in the investigation of the nature of 
the low-lying dipole strength. Recent studies include the application of the Skyrme Hartree- 



Fock + quasiparticle RPA with phonon coupling 



theory 



HQ, 



the time- dependent density-matrix 



jj, the continuum linear response in the coordinate-space Hartree-Fock Bogoliubov 



formalism 



|, the quasiparticle phonon model [2 J3, the relativistic RPA 

and the relativistic quasiparticle RPA (RQRPA) 13, 14]. The PDR is interesting not only 



goliubov 



as a new and exotic mode of nuclear excitations, but it also plays an important role in 
predictions of neutron capture rates in the r-process nucleosynthesis, and consequently in 
the calculated elemental abundance distribution. Even though the El strength of the PDR 
is small compared to the total dipole strength, the occurrence of the PDR significantly 
enhances the radiative neutron capture cross section on neutron-rich nuclei, as shown in 



or the whole nuclear chart 
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it has been shown that 



recent large-scale QRPA calculations of the El stren gth 

In the analysis based on the relativistic (Q)RPA 
in neutron-rich nuclei the electric dipole response is characterized by the fragmentation of 
the strength distribution and its spreading into the low-energy region. In light nuclei, the 
low-lying dipole strength is not collective and originates from non-resonant single-neutron 
excitations. In medium-heavy and heavy neutron-rich nuclei, on the other hand, some of 
the low-lying dipole states display a more distributed structure of the RQRPA amplitudes. 
The study of the corresponding transition densities and velocity distributions revealed the 
dynamics of the neutron pygmy dipole resonance. 

In this work we employ the relativistic QRPA in the study of the evolution of low-lying 
dipole strength in proton-rich nuclei. Because the proton drip-line is much closer to the 
line of /^-stability than the neutron drip-line, bound nuclei with an excess of protons over 
neutrons can be only found in the region of light Z < 20 and medium mass 20 < Z < 50 
elements. For Z > 50, nuclei in the region of the proton drip-line are neutron-deficient rather 
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than proton-rich. In addition, in contrast to the evolution of the neutron skin in neutron-rich 
systems, because of the presence of the Coulomb barrier nuclei close to the proton drip-line 
generally do not exhibit a pronounced proton skin, except for very light elements. Since 
in light nuclei the multipole response is generally less collective, all these effects seem to 
preclude the formation of the pygmy dipole states in nuclei close to the proton drip-line. 
Nevertheless, the present analysis will show that proton pygmy dipole states can develop in 
light and medium mass proton-rich nuclei. 

The relativistic QRPA |13J is formulated in the canonical single-nucleon basis of the 
relativistic Hartree-Bogoliubov (RHB) model and is fully self-consistent. For the interaction 
in the particle-hole channel effective Lagrangians with nonlinear meson self-interactions or 
density-dependent meson-nucleon couplings are used, and pairing correlations are described 
by the pairing part of the finite-range Gogny interaction. Both in the particle-hole and 
pairing channels, the same interactions are used in the RHB equations that determine the 
canonical quasiparticle basis, and in the matrix equations of the RQRPA. This feature is 
essential for the decoupling of the zero-energy mode which corresponds to the spurious 
center-of-mass motion. In addition to configurations built from two-quasiparticle states of 
positive energy, the RQRPA configuration space contains quasiparticle excitations formed 
from the ground-state configurations of fully or partially occupied states of positive energy 
and the empty negative-energy states from the Dirac sea. In the present analysis of the 
dipole response of proton-rich nuclei in the fully self-consistent RHB+RQRPA framework, 



the density-dependent effective interaction DD-ME1 Jj| is employed in the p/i-channel, and 
the finite range Gogny interaction D1S [18] in the pp-channel. 

In Fig.^we display the RQRPA dipole strength distributions in the N=20 isotones: 40 Ca, 
42 Ti, 44 Cr, and 46 Fe. The electric dipole response 



B T (EJ,uo v ) ' 



2Ji + l 

+ (-i)'""-v +J c/° (fi'\\Q T jM k«„v+ (-i)Vv) 



2 



(1) 

is calculated for the isovector dipole operator 

„ N z Z N 

= N + Z ^ rp ^ x,i ~~ N + Z ^ rnYl ^ ' ( 2 ) 

For each RQRPA energy u u , X u and Y u denote the corresponding forward- and backward- 
going two-quasiparticle amplitudes, respectively, v » and Un are the occupation numbers of 



the single-particle levels in the canonical basis. The discrete spectra are averaged with the 
Lorentzian distribution 

*(3=WW0 /) (B _ffi +r , /4 , (3) 

with T = 1 MeV as an arbitrary choice for the width of the Lorentzian. 

The dipole strength distributions in Fig. ^ are dominated by the isovector giant dipole 
resonances (GDR) at « 20 MeV excitation energy. In these relatively light systems the 
GDR still exhibits pronounced fragmentation. With the increase of the number of protons, 
low-lying dipole strength appears in the region below the GDR and, for 44 Cr and 46 Fe, a 
pronounced low-energy peak is found at ~ 10 MeV excitation energy. What is the nature 
of this low-lying dipole state? In the lower panel of Fig. ^ we plot the proton and neutron 
transition densities for the peaks at 10.15 MeV in 44 Cr and 9.44 MeV in 46 Fe, and compare 
them with the transition densities of the GDR state at 18.78 MeV in 46 Fe. Obviously the 
dynamics of the two low-energy peaks is very different from that of the isovector GDR: 
the proton and neutron transition densities are in phase in the nuclear interior and there is 
almost no contribution from the neutrons in the surface region. The low-lying state does 
not belong to statistical El excitations sitting on the tail of the GDR, but represents a 
fundamental mode of excitation: the proton electric pygmy dipole resonance (PDR). 

In Fig. El we analyze the RQRPA structure of the dipole response in 46 Fe. This nucleus is 
located at the proton drip-line (see Ref. for a recent review on the limits of nuclear stabil- 
ity) and, very recently, for the first time evidence for ground-state two-proton radioactivity 
was reported in the decay of 45 Fe 0,|2l|. The discrete unperturbed Hartree-Bogoliubov and 
full RQRPA El spectra of 46 Fe are shown in the left panel. Since the interaction is repulsive 
in the isovector channel, one expects that the inclusion of the residual interaction will result 
in the shift of the full RQRPA spectrum to higher energies with respect to the unperturbed 
spectrum. This is clearly the case for the states in the GDR region. The lowest group of 
states at ~ 10 MeV, however, gains strength and is shifted to lower energy. Obviously the 
nature of these states is not isovector. The QRPA structure of these states is shown in the 
panel on the right, where for the three low-lying states at 9.02, 9.44, and 10.26 MeV, as 
well as for the strongest state in the GDR region at 18.78 MeV, we plot the corresponding 
QRPA amplitude of proton and neutron 2qp configurations 
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^2qp — ^2qp ^2qp i (4) 



with the normalization condition 

£U> = i. ( 5 ) 

For each of the four dipole states, in addition to the excitation energy we have also included 
the corresponding B(E1) value. The QRPA amplitudes are shown in a logarithmic plot as 
functions of the unperturbed energy of the respective 2gp-configurations. Only amplitudes 
which contribute more than 0.01% are shown, and we also differentiate between proton 
and neutron configurations. We note that, rather than a single proton 2qp excitation, the 
low-lying states are characterized by a superposition of a number of 2qp configurations. Ob- 
viously the pygmy states display a degree of collectivity that can be directly compared with 
the QRPA structure of the GDR state at 18.78 MeV. In addition, proton 2qp configurations 
account for « 99% of the QRPA amplitude of the pygmy states, whereas the ratio of the 
proton to neutron contribution to the GDR state is ~ 2. For the GDR states the 2qp con- 
figurations predominantly correspond to excitations from the srf-shell to the /p-shell. The 
structure of the pygmy states, on the other hand, is dominated by transitions from the I/7/2 
proton state at -0.17 MeV, and from the 2p 3 / 2 proton state at 3.60 MeV (this state is only 
bound because of the Coulomb barrier). The energy weighted sum of the strength below 11 
MeV excitation energy corresponds to 2.5% of the Thomas- Reiche-Kuhn sum rule (TRK). 

Another example of particularly pronounced proton PDR are the proton-rich isotopes 
of Ar. In the left panel of Fig. El we display the RHB+RQRPA electric dipole strength 
distribution in 32 Ar. In addition to the rather fragmented GDR structure at ~ 20 MeV, 
prominent proton PDR peaks are obtained just below 10 MeV. For the four states at 8.30, 
8.94, 9.36, and 9.60 MeV, which constitute the pygmy structure, in Fig. 0] we plot the 
proton and neutron transition densities, and the distributions of the RQRPA amplitudes, 
in comparison with the GDR state at 18.01 MeV. In contrast to the isovector GDR, the 
proton and neutron transition densities of the pygmy states are in phase. The RQRPA 
amplitudes of the low-lying states present superpositions of many proton 2qp configurations, 
with the neutron contributions at the level of 1%. The dominant configurations correspond 
to transitions from the proton states 1^3/2 at -2.09 MeV and 2sy2 & t -4.07 MeV. The 
corresponding B(E1) values are included in Fig.HJ The energy weighted sum of the strength 
below 10 MeV takes 5.6% of the TRK sum rule. In the right panel of Fig. El we display the 
mass dependence of the centroid energy of the pygmy peaks and the corresponding values 
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of the integrated B(E1) strength below 10 MeV excitation energy. In contrast to the case 
of medium-heavy and heavy neutron-rich isotopes, in which both the PDR and GDR are 
lowered in energy with the increase of the neutron number 2], i n proton-rich isotopes the 
mass dependence of the PDR excitation energy and B(E1) strength is opposite to that of 
the GDR. The proton PDR decreases in energy with the growth of the proton excess. This 
mass dependence is intuitively expected since, as we have shown, the proton PDR mode is 
dominated by transitions from the weakly-bound proton orbitals. As the proton drip-line 
is approached, either by increasing the number of protons or by decreasing the number 
of neutrons, due to the weaker binding of higher proton orbitals one expects more inert 
oscillations, i.e. lower excitation energies. The number of 2qp configurations which include 
weakly-bound proton orbitals increases towards the drip-line, resulting in an enhancement 
of the low- lying B(E1) strength. 

In conclusion, we have employed the self-consistent RHB model and the RQRPA in the 
analysis of the low-energy dipole response in proton-rich nuclei close to the proton drip-line. 
For nuclei with more protons than neutrons, in the chain of N=20 isotones the RHB+RQRPA 
calculation predicts the occurrence of an exotic mode in the low-energy region below 10 MeV: 
the proton electric pygmy dipole resonance. The analysis of the corresponding proton and 
neutron transition densities shows that proton PDR states correspond to the oscillation of 
the proton excess against an approximately isospin-saturated core. The RQRPA amplitudes 
of the pygmy states display a superposition of many proton 2qp configurations, whereas 
the contribution of neutron excitations is at the level of 1%. The dominant configurations 
include weakly-bound proton orbitals, and also the states which are bound only because 
of the presence of the Coulomb barrier. In contrast to the rather fragmented GDR, the 
PDR strength is concentrated in a narrow region of excitation energies well below the region 
of giant resonances, and the energy weighted sum of the PDR strength corresponds to a 
few percent of the TRK sum rule. With the decrease of the number of neutrons, for the 
chain of proton-rich Ar isotopes the proton PDR is lowered in energy, and the integrated 
B(E1) strength in the low-energy region increases accordingly. For heavier nuclei the proton 
drip-line is located in the region of neutron-deficient, rather than proton-rich nuclei, and 
therefore no low-lying dipole strength is calculated in medium-heavy and heavy nuclei close 
to the proton drip-line. 
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FIG. 1: The RHB+RQRPA isovector dipole strength distributions in the N=20 isotones, calculated 
with the DD-ME1 effective interaction. For 44 Cr and 46 Fe the proton and neutron transition 
densities for the main peak in the low-energy region are displayed in the lower panel and, for 46 Fe, 
the transition densities for the main GDR peak. 
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FIG. 2: The discrete Hartree-Bogoliubov and full RQRPA El spectra in 46 Fe (left panel). For 

the main peaks in the low-energy region and in the region of the isovector GDR, the distributions 
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of the RQRPA amplitudes are shown as functions of the unperturbed energy of the respective 




FIG. 3: The RHB+RQRPA isovector dipole strength distribution in 32 Ar (left panel). For the argon 
isotopes, the mass dependence of the centroid energy of the pygmy peak and the corresponding 
values of the integrated B(E1) strength below 10 MeV excitation energy are shown in the panel on 
the right. 
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FIG. 4: Transition densities and RQRPA amplitudes for selected 1~ states in 32 Ar. For the main 
peaks in the low-energy region and in the region of the isovector GDR, the proton and neutron 
transition densities are shown in the left panel. The corresponding distributions of the RQRPA 



